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ABSTBUCT 


Ihe  lEH  Ccapany*s  Continasus  Systeas  Hodeling  Program 
was  used  to  siaulata  th«  longitudiaal  flight  control  sysrea 
of  the  P/&-18  aircraft.  The  aodel  is  intended  for  ase  in 
investigations  cf  aircraft  response  to  flight  conditions 
which  approach  spin  or  stall  and  is  restricted  to  the 
aatoaatic  flaps  ap  (AFO)  flight  node.  Prograa  outputs 
include  stabllator  deflection,  leading  and  trailing  edge 
flap  positions,  and  cross-axis  interconnect  signals. 
Various  stick  forces,  aoticn  sensor  Inputs,  and  air  pressure 
inputs  were  sisulated  to  produce  transient  control  surface 
responses.  These  coaputcr  generated  responses  exhibited 
characteristics  corresponding  to  predicted  aircraft  control 


surface  mowescnts 
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I.  IHTRODaCTIOM 


The  navy  has  experienced  the  loss  of  nuneroas  aircraft 
during  recent  years  due  to  unintentional  departure  from 
controlled  flight.  The  increased  cost  and  complexity  of 
modern  aircraft  utilizing  fly-by-wire  flight  control  systems 
have  placed  a  renewed  emphasis  on  understanding  the 
performance  characteristics  of  these  aircraft,  especially 
naar  the  limits  cf  the  flight  envelope.  The  purpose  of  this 
thesis  is  to  investigate  a  method  of  non-real  time  computer 
simulation  of  the  longitudinal  flight  control  system  of  the 
F/&-18  aircraft.  Future  thesis  researchers  at  the  Naval 
Postgraduate  School  will  complete  corresponding  simulations 
of  the  lateral  and  directional  control  systems.  The 
combination  of  these  simulations  with  an  existing 
aerodynamic  simulation  program  will  yield  a  complete 
aircraft  stability  and  control  model.  The  primary  purpose  of 
the  model  is  to  investigate  methods  of  designing  control 
augmentation  systems  which  actively  inhibit  or  prevent 
departures  from  controlled  flight.  Othar  uses  of  this  model 
would  include  the  capabilities  to  test  new  programmable 
memory  configurations,  to  evaluate  new  components  such  as 
ootimal  observers,  to  sirulate  degraded  flight  conditions 
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(such  as  a  daoaged  flight  control  surface)  and  to  recreate 
flight  conditions  daring  post  accident  investigations. 

The  sethod  which  was  chosen  for  acconplishment  cf  the 
objectives  of  this  -*'hesis  was  the  Continuous  Systeas 
Modeling  Prograa  (CSHP)  ,  developed  by  the  IBM 
Coapany  [Bef.  1].  CSHP  is  a  software  package  designed  to 
siaulate  djnaaic  systess  described  in  teras  of  differential 
equations  and  block  diagrams  noraally  encountered  in  systeas 
theory.  CSHP  allows  prograaaing  flexibility  through  the  use 
of  thirty-four  pre-prograaaed  functional  blocks  which  are 
siailar  tc  FOBTBAN  subroutines.  These  blocks  provide  rapid 
access  to  aatheaatical  functions,  switching  functions, 
signal  sources,  logic  functions  and  FOBTSAN  functions,  since 
this  thesis  represented  the  first  atteapt  at  the  Naval 
Postgraduate  School  to  accurately  aodel  the  flight  control 
systea  of  a  modern,  highly  augaented  tactical  aircraft  it 
was  deeaed  important  to  concentrate  on  the  physical  systeas 
rather  than  becoae  Involved  in  the  coaplezities  of  numerical 
analysis. 

Alternatives  to  CSHP  which  were  considered  included 
analcg  prograaaing  and  FCBTRAN  prograaaing.  Analog 
prograaaing  was  not  selected  because  it  is  less 
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rapresentative  cf  1:he  systems  to  be  modeled  and  it  is  less 
accurate  than  CSHP.  The  concept  of  programming  directly  in 
PORTBAN  uas  carefully  considered,  since  FORTRAN  is  the 
source  language  for  CStIP  the  capabilities  of  cSHP  are  a 
subset  of  the  capabilities  of  FORTRAN  itself.  Additionally, 
CSHP  has  resxrictions  on  the  number  of  allowable  statements, 
constants,  variables  and  ether  parameters.  Dnlike  FORTRAN, 
when  functional  blocks  in  CSHP  axe  used,  the  programmer  has 
no  direct  control  of  mathematical  operations  internal  to  the 
functions.  The  primary  reasons  for  which  CSHP  was  selected 
were  its  simplified  input  statements,  output  staxeaents  and 
program  control  statements  which  facilitated  rapid  program 
writing  and  testing.  Additionally,  the  automatic  time  and 
amplitude  scaling,  data  formatting,  and  compatibility  with 
graphic  display  devices  which  CSHP  provides  are  well  suited 
for  prototype  pregram  development. 


II 


k.  PBOGRAH  0B6AM1ZATI0M  AHD  NON ENCLATOBE 

A  detailed  description  of  the  flight  control  laws  of  the 
F/A-18  aircraft  is  given  in  the  McDonnell  Aircraft  Company's 
system  design  report  [Bef.  21.  The  comtater  program 
developed  in  this  thesis  is  based  upon  figure  16.  1  of  this 
report,  entitled  P/A- 18  Lcngitudinal,  Mechanical,  CAS,  and 
DEL  Control  Lav  Mechanization.  This  figure  connains  six 
pages  of  block  diagrams  depicting  generation  of  longitudinal 
control  signals  which  are  valid  for  all  aircraft 
configurations  and  failure  modes.  The  version  of  flight 
control  program  incorporated  in  the  flight  control  computer 
programmable  read  only  memory  (PROM)  utilized  in  this 
simulation  is  8.2.1,  dated  August  31,  1982. 

A  brief  discussion  of  the  P/A- 18  flight  control  system 
is  necessary  to  facilitate  a  discussion  of  control  law 
modeling.  All  computations  of  control  laws  are  accomplished 
independently  by  four  channels  of  digital  computation. 
Primary  ccntrol  in  the  pitch  axis  is  provided  by  symmetric 
deflection  cf  the  horizontal  stabilators.  Pull  span  leading 
edge  flaps  and  trailing  edge  flaps  are  scheduled  to  provide 
maximum  lift  tc  drag  ratio  during  maneuvering,  high  angle 
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of  attack  and  cruise  configurations.  Roll  control  is 
accoaplished  by  ailerons,  differential  stabilators  and 
differential  leading  and  trailing  edge  flap  deflection. 
Directional  control  is  eaintained  by  dual  rudders.  The 
thesis  coeputei  program  which  is  contained  in  Appendix  A 
simulates  the  output  of  one  channel  of  digital  computation 
and  calculates  the  angular  positions  of  the  stabilators, 
leading  edge  flaps  and  trailing  edge  flaps.  This  simulation 
does  not  calculate  rudder  or  aileron  deflection,  however, 
all  electrical  signals  required  for  cross  axis  interconnects 
are  provided. 

The  task  cf  programming  the  information  given  in  the 
longitudinal  ccntrol  law  mechanization  schematic  was 
simplified  by  two  means.  First,  the  program  restrictions  and 
assumptions  to  conditions  cf  flight  which  are  discussed  in 
part  B  were  applied.  As  a  result  sections  of  figure  16.1 
which  apply  tc  mechanical  ccntrol  laws  and  spin  modes,  for 
example,  were  deleted.  This  reduced  the  number  of  schematic 
blocks  to  be  modeled  by  approximately  one  third.  The  second 
simplification  arose  through  a  system  of  nomenclature  in 
which  nine  control  paths  were  defined  in  order  to  limit  the 
number  of  input  and  output  signals  for  any  specific  path. 


When  ccnbined,  these  nine  paths  form  the  total  longitudinal 
control  lav  aechanization . 

Figure  2.1  is  the  overall  longitudinal  control  signal 
block  diagras  vblch  was  used  in  this  siaulation.  It  was 
derived  by  applying  all  program  assumptions  or  restrictions 
outlined  in  part  B  to  figure  16.1  of  the  system  design 
re  port. 

Block  diagrams  depicting  the  logic  development  of  the 
component  paths  are  included  in  Appendix  A.  The  nomenclature 
for  each  control  path  which  is  given  in  Table  I  is  peculiar 
to  this  simulation  program. 

Table  II  lists  nomenclature  for  control  signal  groups  which 
are  common  to  both  the  BcDonnell  schematic  and  to  this 
siaulation.  Signals  with  common  prefixes  are  numbered 
consecutively  in  the  feed  forward  direction.  The  primary 
feed  forward  input  to  the  CSHP  siaulation  is  pilot  stick 
force.  Feedback  signals  include  Ditch  rate,  roll  rater  yaw 
rater  angle  of  attack,  ncraal  acceleration  and  differential 
control  surface  commands. 

B.  FBOGBAB  BE5TB1CTI0NS  AHD  ASSOHPTIONS 

The  task  of  obtaining  a  working  computer  program  to  meet 
the  thesis  objectives  did  not  require  the  siaulation  of  all 
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possible  conditicns  of  flight.  Thus*  an  assuaed  aircraft 
configuration  and  flight  condition  led  to  aany 
siaplifications  in  the  acdel.  Each  aajor  assuaption  is 


discussed  below.  Should  future  researchars  desire  to 


construct  a  acre  general  acdel*  additional  prograa  logic 
paths  could  be  included  in  this  CSHP  siaularion  without 
requiring  sajor  rewisions  to  the  prograa. 
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Issttapticss: 

1.  The  aircraft  flight  control  electronics  set  (FCSS)  is 
operating  noreally.  Friaary  coapcnents  of  this  set  include 
flight  control  ccaputers*  pitch*  roll  and  yaw  rate  gyros* 
acceleroaeter  asseablies  and  flight  control  panels.  The  FCES 
contains  logic  seguences  for  failure  defection  and 
corrective  control  lav  iipleaentation  vhich  were  oaaitted 
froa  the  thesis  coapurer  prograa , 

2.  The  aircraft  is  under  inner  loop  control.  In  rhis 
aode  pilot  stick  force  foras  the  priaar/  control  input. 
Autopilot  functions  such  as  heading  hold*  roll  or  pitch 
attitude  hold*  and  altitude  hold  are  not  in  operation. 

3.  The  pitch  control  augaentation  systea  (CAS)  is 
oroviding  longitudinal  control.  Backup  control  systeas  vhich 
are  not  in  operation  include  the  aechanical  backup  aode  and 
the  direct  electric  link  (CEL)  aode  which  provides  an  open 
loop  signal  frca  the  pilct  stick  position  sensor  to  the 
stabilator  servoactuator . 

4.  The  aircraft  is  operating  in  the  u?  and  away  flight 
phase  in  the  auto  flaps  up  (AFO)  configuration.  This  phase 
requires  that  the  flap  switch  in  the  cockpit  be  in  the  AOTO 
position  or  that  the  calibrated  airspeed  be  greater  than  24  3 
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knots.  In  ths  AFU  node  gain  schedules  are  optimized  for 
combat  maneuvering  characteristics  in  the  low  to  aid  dynamic 
pressure  regime.  Addit icnally,  a  trim  integrator  controls 
load  factor  and,  at  angles  of  attack  above  22  degrees, 
proportional  nose  down  commands  are  introduced. 

5.  The  gain  switch  is  in  the  normal  position.  The 
override  position  of  this  switch  causes  flight  computers  to 
use  fixed  values  for  air  data  schedules  and  a  predetermined 
angle  of  attack  in  control  lav  computations.  The  normal 
oosition  of  the  gain  switch  allows  measured  values  of  air 
data  and  angle  of  attack  to  be  used. 

6.  The  aircraft  is  operating  with  weight  off  the  wheels, 
soeed  brake  in,  and  with  no  external  stores. 

7.  Anti  spin  functions  are  not  simulated. 
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A.  AIB  DATA  SCBEDOLES 


Air  data  schednlas  are  functions  of  static  pressure 
(PS)  «  cospressible  dynaiic  pressure  (QC) «  and  other 
paraaeters  such  as  noraal  acceleration,  angle  of  attack  or 
condition  of  external  stores.  The  longitudinal  control  laws 
contain  inputs  froa  13  different  air  data  schedules  which 
are  listed  in  Table  III. 


TABLE  III 


Air  Data  Schedules 


POBCTiOH 

P12 

P22 

123 

P24 

P27 

P29a 

P|2A 

P37 

140 

P68 

f107 


PiSCBIPlIOH  (xmis) 

Pitch  forward  loop  integrator  gain  schedule  (BI,PS) 
Pader  on  supersonic  coapensatlpn  (QC) 

Stall  aargin  on  pitch  forward  loop  integrator  (AOA) 
Trailing  edge  flap  schedule  (AOA. Rif 
Trailing  edge  flap  schedule  roc  liait) 

LeaMng  edge  flap  schedule  (AOA.fl) 


Leading  edge  flap  schedule  (AOA.BI) 

Leading  edge  flap  schedule  (QC  liait) 

Leading  edge  flap  schedule  (Si  liait) 
Longitudinal  forward  loop  gain  schedule  (QC) 
HZ  liait  on  AOA  feedback  (HZA) 

Pitch  rate  feedback  gain  schedule  (PS.QC,SI) 
Pitch  rate  feedback  gain  schedule  (QC) 
Longitudinal  Inertial  gain  schedule  (QC) 


Each  schedule  controls  an  output  gain  for  a  particular 
purpose.  Function  32A,  for  ezaaple,  yields  a  unifora  initial 
Ditch  acceleration  in  response  tc  sharp  inputs,  with  gain 
decreasing  at  high  walues  of  coapressible  dynaaic  pressure. 
Other  functions  such  as  function  290  are  used  to  deteraine 
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an  upper  or  lover  liait  to  an  input  signal,  siaulation  of 
these  functions  in  the  CSHI  foreat  required  that  Individual 
cosputer  prograas  be  written  to  test  each  air  data  schedule. 
In  all  cases  the  final  result  was  given  in  the  flight 
control  systes  design  report,  usually  in  graphical  form. 
Conversion  to  CSHP  format  was  done  by  extraction  of  data 
points  from  graphs,  use  of  logic  flow  charts  when  available 
and  direct  employment  of  those  mathematical  formulae  which 
were  given.  Cnee  programmed  in  CSHP,  the  tabular  and 
graphical  output  data  was  compared  to  data  given  in  the 
design  report  for  the  same  test  parameters.  Thus,  the  output 
of  each  function  was  independeirtly  verified  before  the 
function  was  Included  in  the  longitudinal  control  law 
siaulation. 

The  siaulaticn  is  written  tc  take  advantage  of  CSHP*s 
SORT  and  HGSOBT  capabilities  where  appropriate.  The  NOSOBT 
'option  of  CSHP  is  used  for  conditional  logic  and  branching. 
This  statement  allows  the  user  flexibility  in  creating 
sections  cf  the  program  in  which  ordinary  FORTBhN  rules  can 
be  used.  In  HOSORT  sections  intermediate  variables  were 
defined,  comparisons  made  and  appropriate  branching  was 
executed.  The  program  statements  were  returned  to  SORT 
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foraat  as  soon  as  conditional  logic  was  no  longer  required. 
kn  early  prcblea  revealed  that  if  two  or  more  SORT  secticns 
are  separated  hy  a  NOSORI  section  inforaacion  will  not  be 


passed  between  the  separated  SORT  sections.  To  prevent 
errors  resulting  froa  this  restriction  froa  occuring  the 
follcwing  decisions  regarding  the  order  of  program 
statements  were  made. 

1.  The  nuater  of  NOSOBI  secticns  should  be  ainiaized. 

2.  NOSOBT  sections  should  be  located  close 
together , allowing  fewer  and  larger  SORT  sections. 

3.  Sacros  should  be  utilized  when  possible.  Macros , 
which  are  similar  to  subroutines,  ar'e  discussed  in  part  0 
under  frequency  averagers. 

The  task  of  arranging  prograa  stateaents  to  ainiaize  the 
number  of  BOSOBT  sections  was  greatly  siaplified  through  the 
use  cf  the  "output  variable  sequence"  tabularion  which  is 
produced  as  part  of  the  CSBP  standard  output.  The  final 
nuaber  of  HOSOBT  stateaents  in  the  computer  simulation  could 
have  been  further  reduced  by  this  aethod,  however,  this 
would  have  required  the  acveaent  of  large  program  blocks  and 
caused  a  deviation  froa  the  logic  oath  used  in  prograa 
development.  For  exaaple,  the  NOSOBT  sections  of  functions 
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12  aad  40  coaid  have  been  conbined.  This  would  hsive 
Increased  difficulty  in  progran  debugging  and  nade  areas  of 
the  prograa  which  require  coabinational  logic  less  apparent 
to  readers. 

Appendix  C  contains  docuaentation  for  the  air  data 
schedules.  Bach  schedule  is  included  in  a  complete  computer 
program  which  produces  tabular  and  graphical  data  to  match 
the  characteristics  of  figures  in  the  design  report.  In 
incorporating  these  air  data  schedules  into  the  longitudinal 
control  lav  simulation  all  schedules  except  for  function  24 
were  placed  before  the  main  computational  body  of  the 
program,  function  24,  the  trailing  edge  flap  schedule, 
requires  conditional  logic  and  contains  computed  angle  of 
attack  as  an  input.  Since  computed  angle  of  attack  is 
generated  in  the  normal  accelerometer  path,  function  24  was 
placed  immediately  before  the  trailing  edge  flap  path 
computations. 

B.  PADBRS,  ABAIOS  TO  DIGITAL  CONVERTERS,  BECHAHICAL 

BEBAROOTS 

The  purpcse  of  faders  is  to  eliminate  large 
discontinuities  in  gain,  permitting  gradual  change  in  output 
from  old  values  to  new  new  values  at  a  desired  sample  rate. 
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Pad«zc  arc  located  in  Figure  16.1  at  the  outputs  of  air  data 
schedoles  and  gain  schedules  which  are  dependent  upon 
aircraft  configuration.  Discontinuities  aay  arise  as  a 
result  of  a  change  in  physical  seasurenents  such  as  dynamic 
pressurOr  change  in  aircraft  configuration  such  as  speed 
brake  extension^  or  change  in  mode  of  flight.  The  lower 
limit  for  signal  HP4,  for  example,  changes  from  a  gain  of 
>50.0  to  0.0  as  a  result  cf  spin  entry.  Since  the  thesis 
computer  program  assumes  that  electrical  signals  vary 
smoothly  and  is  restricted  to  up  and  away  controlled  flight, 
faders  are  modeled  in  freguency  only. 

inalog  to  digital  conversion  and  frequency  matching  are 
obtained  by  a  sample  and  held  process.  Sampling  tines  are 
generated  by  the  CSdP  functional  block  IHPOLS  which  produces 
a  time  series  of  unit  impulse  functions  with  a  specified 
start  time  and  period.  Since  these  pulse  trains  are  used  in 
several  areas  cf  the  program,  impulse  functions  of  20,  40 

and  80  hertz  are  included  immediately  following  the  air  data 
schedules.  The  zero  order  hold  function  ZHOLO  keeps  signal 
gain  constant  throughout  the  pulse  period. 

1  mechanical  breakout  force  of  two  pounds  is  modeled  in 
the  pilot  stick  input  path.  Vhen  large  stick  forces  are 
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aopli«d,  the  resultant  signal  is  calculated  with  a  redaction 
in  magnitude  c£  two  pounds.  This  mcdelinq  is  accomplished  in 
eSHP  by  the  deadspace  (OEADSP)  functional  block.  Figure  16.1 
of  the  design  report  depicts  electrical  signal  "deadbands'* 
which  are  ccnceptually  identical  to  mechanical  breakouts. 
Appendix  D  contains  a  program  which  displays  deadspace 
outputs  in  tabular  and  graphical  form. 


C.  ALIASING  FILTERS  AND  SIGNAL  LIMITERS 

Ihe  longitudinal  control  laws  cf  the  F/A-18  include  five 
aliasing  filters  which  are  modeled  as  first  and  second  order 
Laplace  transforms.  The  first  order  transform  is  of  the 
type  A/(Bsel)  and  is  in  tie  pilot  stick  input  path.  This  lag 
type  filter  with  one  real  pole  is  converted  tc  CSMP  format 
by  the  functional  block  REALPL.  A  required  initial  condition 
is  the  value  cf  the  output  signal  when  time  is  zero.  This 
may  be  determined  arbitrarily  by  the  user,  but  was  set  to 
zero  for  this  simulation. 

Second  order  filters  are  present  in  the  pilot  stick 
input  path#  pitch  rate  gyro  path#  angle  of  attack  sensor 
oath  and  normal  accelerometer  path.  Each  is  of  the  form 
A/ (Bs^ i-Cs^O)  and  represents  an  underdamoed  system.  The  CSMP 


functional  block  for  complex  poles  (CMPXPL)  is  used. 
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Initial  conditions  are  tie  value  of  the  output  signal  and 
time  rate  of  change  of  tha  output  signal  when  time  equals 
zero.  Natural  frequencies  varied  from  4.34  hertz  in  the 
pilot  stick  input  path  to  33.3  hertz  in  the  angle  of  attack 
sensor  path.  The  appendix  contains  a  computer  program  which 
demonstrates  tfce  CSMP  outputs  of  both  first  and  second  order 
aliasing  filters.  Unit  step  inputs  were  introduced  to  each 
filter  to  generate  transient  responses.  The  first  order 
filter  produced  an  exponential  rise  to  steady  state  with  the 
correct  tine  constant.  Characteristics  of  the  second  order 
filters  such  as  rise  tine,  peak  time,  maximum  overshoot,  and 
settling  tine  compared  favorably  with  theoretical 
results  CB«f.  3]. 

Signal  limiters  restrict  the  maximum  or  minimum  values 
of  an  output  signal.  Stabilator  surface  deflection,  for 
example,  is  limited  to  10.5  degrees  trailing  edge  down  and 
24  degrees  trailing  edge  up.  The  CSMP  functional  block  LIMIT 
allows  direct  specification  of  lower  and  upper  signal 
limits. 

D.  FBEQOEHCT  JIVEBAGEBS  AND  BATE  LIHITER5 

A  characteristic  of  the  F/A- 18  flight  control  system  is 
that  various  signal  paths  operate  at  frequencies  of  29,  40 


or  80  hertz.  Nhen  signals  are  combined  mathematically,  the 
inputs  are  first  converted  to  a  ccmmon  frequency.  Normal 
accelerometer  path  signals  are  computed  at  40  hertz  and 
combined  with  outputs  of  the  forward  integrator  path,  which 
operates  at  20  hertz.  In  this  case  a  20  to  40  hertz  averager 
is  present  between  the  integrator  path  and  the  summing 
junction  connected  to  the  normal  accelerometer  path.  The 
algorithm  used  for  the  20  to  40  hertz  averager  is  based  on  a 
orocedure  given  in  the  Flight  Control  Electronic  System 
Report  [Ref.  4].  The  averager  was  required  to  generate 
signal  values  at  twice  tie  rate  of  the  incoming  pulses  by 
linear  interpclation  between  amplitudes  of  the  two  previous 
signals  at  20  hertz.  The  formulas  used  to  generate 
intermediate  signals  were: 

Z40  =  Z40Z 1  ♦  DEL  (3.1) 

DEI  =  (Z20  -Z20Z1)  /  2.0  (3.2) 

where : 

Z4Q  s  current  value  of  the  40  hertz  signal 
Z40Z1  *  previous  value  of  the  40  hertz  signal 
Z20  »  current  value  of  the  20  hertz  signal 

Z20Z1  a  previous  value  of  the  20  hertz  signal. 


Conditional  logic  was  usad  to  keep  the  output  signal 
equal  to  the  input  signal  at  times  when  tha  20  hertz  impuse 
function  was  equal  to  one.  The  algcritho  was  initialized  by 
letting  "previous  values"  be  equal  to  input  values  when  time 
equals  zero. 

Ihe  presence  of  numerous  frequency  averagers  in  the 
control  laws  would  have  required  the  repetition  of  many 
statement  blocks  without  the  use  of  CSMP  program  MACROS, 
which  are  similar  tc  FCETBAM  subroutines.  Frequency 
averagers  in  tbe  thesis  computer  simulation  were  included  in 
MACROS  and  placed  at  the  beginning  of  the  program.  A  MACRO 
may  be  used  several  times  within  a  program.  Input  and  output 
variables  are  given  dummy  names,  yet  the  MACRO  is  invoked 
with  a  unique  name 'which  is  assigned  in  a  function 
definition  statement.  A  limitation  to  the  use  of  MACROS  is 
that  variables  which  are  defined  in  MACRO  structure 
statements  are  not  available  for  output  unless  they  are 
designated  as  arguments  in  the  function  definition 
statement.  Additionally,  certain  functional  blocks  such  as 
BEALPL,  CHPXPL  and  IMT6RL  cannot  be  used  as  arguments  or  as 
parts  of  a  MACRO  structure  statement. 
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The  PBOCEDORE  function  of  CS HP  was  used  in  each  HACRO  to 
cause  stateaents  to  be  executed  in  the  order  of  their 
appearance.  Each  structure  statement  of  a  MACRO  will  be 
individually  sorted  unless  PBOCEDORE  is  specified.  All 
statements  contained  within  a  PBOCEDORE  function  are  treated 
as  a  single  block  which  can  be  moved  but  not  rearranged  by 
the  CSHP  translator. 

FORTBAR  subprograms  were  not  used  in  this  simulation 
because  the  size  of  the  Icngitudinal  portion  of  the  program 
did  not  approach  the  maziaum  limits  of  CSHP.  Inclusion  of 
the  lateral  and  direc'^'ional  control  systems  in  this  program 
will  require  that  measures  be  taken  to  remain  within  the 
allowable  number  of  structure  statements,  NOSORT  sections 
and  HACROS.  The  CSHP  translator  is  not  used  to  process 
FORTRAN  subprogram  stateients.  Since  the  number  of 
subprogram  statements  is  net  counted  the  overall  size  of  the 
program  may  be  increased.  The  capabilities  of  the  computer 
system  library  may  be  utilized  by  subprograms  through  the 
use  of  the  CALL  statement.  The  CALL  statement  must,  however, 
be  included  in  a  NOSORT  or  PBOCEDORE  section.  A  method  to 
invoke  subprograms  exists  which  does  not  require  NOSORT  or 
PBOCEDORE  sections,  but  it  is  valid  only  for  two  or  more 
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output  variables  and  it  requires  a  specific  format  for 
arqusents.  A  final  restriction  to  the  use  cf  FOBIBAM 
subpregrass  is  that  certain  CSHP  functional  bloclcs  such  as 
ZHOLO,  IHPOIS,  and  CHPXPL  are  not  allowable. 

The  longitudinal  control  laws  contain  rate  limiters 
which  operate  at  frequencies  of  20  hertz  and  80  hertz  to 
restrict  the  speed  of  leading  edge  flap  movement.  The 
algorithm  for  these  limiters  compares  the  value  of  each 
inccmir.g  signal  with  the  value  of  the  previous  output 
signal.  The  magnitude  of  the  difference  between  these 
signals  is  processed  by  the  LIN  IT  functional  block  which 
generates  the  current  output  signal.  The  20  hertz  rate 
limiter,  for  example,  allcws  a  *  naziaum  change  in  output 
signal  value  cf  0.9  degrees  during  each  period  of  0.05 
seconds.  A  listing  of  the  program  which  was  used  to  test  the 
rate  limiter  is  contained  in  appendix  D.  Correct  outputs 
were  observed  fer  both  increasing  and  decreasing  input 
signals. 

B.  CI6ITA1  TILTEBS,  DIGITAL  TO  ANALOG  CONVEBTEBS, 

SEBVOEECBANISHS 

Three  types  of  digital  filters  are  used  in  the 
longitudinal  ccntrol  lavs.  Z  filter  number  P2  is  a  lead‘‘lag 


type  controller  in  the  pitch  rate  gyro  path  which  operates 
at  20  hertz.  Z  filter  number  ?9,  the  forward  loop 
integrator,  operates  at  20  hertz  and  compares  the  aircraft 
response  to  the  maneuwer  command.  The  output  signal  driwes 
the  stabilator  servoactuatcr  to  reduce  the  maneurer  error  to 
zero.  This  allows  the  aircraft  to  to  be  automatically  kept 
in  a  hands  off  condition  since  the  forward  integrator 
eliminates  onccmvanded  normal  acceleration.  Z  filter  number 
P8  is  a  structural  notch  filter  which  operates  at  80  hertz. 
It  attenuates  aeroelastic  bending  which  is  detected  by  the 
moticn  sensors. 

In  the  simulation  each  filter  was  developed  in  its  most 
aeneral  form  for  inclusion  in  a  MACRO.  The  following 
equation  for  lead-lag  filter  number  P2  is  given  in  the 
Schematic  Cesign  Report  C  5]. 

PS^  (UgR11*(1«gKl2))Z  -  (PK1U1)»(1-?K12) 

PR3  ”  Z  -  (1-PK12)  (3.3) 

It  is  mcdeled  in  the  thesis  simulation  as 

FOOT  AZ  •  B 

FIH  Z  -  C  (3.4) 

where  A,  B  and  C  are  constants.  The  right  shifting  and 
linearity  properties  of  the  z  transform  are  used  to  solve 

explicitly  for  the  variable  FOOT  TBef.  6]. 

30 


POUT 

Pin 


A  -  B  (2-  1) 
1.0  -  C(Z-1) 


(3.5) 


Cross  aoltiplica':ion  and  rearranges en-t  yields: 

POOT*  A*PI1I  -  B*PIH(Z-1)  ♦  C*POOT(Z-1)  (3.6) 

which  is  described  in  the  sisulation  as: 

POOT*  A*PIH  -  BepINZI  ♦  CePOUTZI  (3.7) 

This  sethod,  which  is  termed  direct  realization 
programming,  was  also  used  in  the  development  of  the  notch 
filter  and  the  forward  loop  integrator.  The  equivalent 
Laplace  transform  for  each  longitudinal  flight  control 
filter  is  listed  in  the  design  report.  This  permitted  a 
cress  check  of  z  filter  performance  which  is  included  in 
aopendiz  D. 

A  specific  method  of  integration  may  be  specified  in  the 
terminal  portion  of  a  CSHP  program.  In  the  case  of  flight 
control  simulation  the  Bunge^Kutta  Fixed  Step  size  (RKSFX) 
method  was  utilized  to  ensure  that  integrations  would  only 
occur  at  the  desired  sampling  rate.  The  highest  sampling 
frequency  in  any  axis  of  the  P/A -18  flight  control  system  is 
80  hertz,  thus  the  CSHP  integration  interval  DELT  was 


specified  as  0.0125  seconds 


The  ▼ariatl€  KEEP  is  used  in  CSHP  to  indicate  that  the 
end  cf  a  valid  integration  step  has  been  reached.  KEEP  is 
set  equal  to  one  when  this  conditicn  is  set.  During  trial  or 
interaediate  integration  steps  KEEP  will  equal  zero.  Each 
HKCBC  contains  conditional  logic  which  allows  calculations 
to  be  perforaed  only  when  KEEP  equals  one. 

Conversion  of  signals  froa  analog  to  digital  fora  in  the 
P/A-18  occurs  as  the  signal  reaches  the  servomechanism.  The 
quantizer  functional  block  (QNTZ  B)  is  employed  to  accomplish 
the  analog  to  digital  siaulation.  The  transfer  functions 
used  by  the  stabilator  and  flap  servomechanisms  are  not 
published  by  the  manufacturer  of  the  aircraft.  In  order  for 
the  thesis  coaputer  program  to  generate  control  surface 
positions  the  response  characteristics  of  a  Parker  Hannifin 
fly-  ty-wire  actuation  system  were  incorporated.  The  selected 
actuators  were  designed  for  use  with  all-digital  flight 
controls  and  are  modeled  as  second  order  systems.  The 
stabilator  transfer  functicn  is 

X  1600 _ 

eI  *  S6S  ♦  1600  (3.8) 

where  Z  is  the  actuator  position  in  degrees  and  Ei  is  the 
position  ccaaand.  The  transfer  function  for  both  leading 
and  trailing  edge  flaps  is 
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X 

e1 


400 

S*+  28S  *4  00 


(3.9) 


The  advantages  of  using  the  second  ordar  nodel  instead  of  a 
first  order  eodel  are  that  the  faster  rise  time  aore  closely 
represents  the  physical  actuators  and  that  the  natural 
frequency  and  daaping  ratio  aay  be  independently  modified. 

P.  FBOGRAH  TESTIMG  METHODS  AND  RESULTS 

The  thesis  computer  program  was  tasted  on  three  levels. 
The  lowest  level  involved  evaluation  of  the  individual 
signal  blocks  cf  figure  16.1  in  the  system  design  report. 
Sections  A  through  B  of  this  chapter  describe  signal  block 
modeling  techniques  and  appendicies  C  through  6  contain 
testing  programs.  To  obtain  verification  of  proper  program 
operation  it  was  desired  to  create  input  signals  of 
realistic  value  which  would  produce  a  time  varying  output. 
The  rate  limiters,  for  example,  were  tested  with  an  input 
signal  which  rose  exponentially  to  a  limiting  value,  then 
decreased  expcrentially.  Output  signals  at  the  desired 
frequency  were  observed  for  incoming  signals  of  positive  or 
negative  slope  and  for  incoming  signals  which  were  within  or 
beyond  the  rate  limit. 
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The  second  level  of  progtaa  testing  involved  the  nine 
signal  paths  which  are  listed  in  Table  1.  In  most  cases  a 
step  function  was  used  as  the  input  siqnal.  Intermediate 
signals  were  observed  to  determine  continuity, 
freguency ,time  constants,  and  conformance  with  limiting 
values. 

The  highest  level  of  program  testing  required  the 
aeneration  of  ccntrol  surface  deflections  for  specific 
combinations  of  pilot  stick  force  and  motion  sensor  feedback 
signals.  Neither  aircraft  flight  data  nor  the  McDonnell 
flight  control  simulation  data  were  available  for  direct 
co'mparison  with  the  outputs  of  the  thesis  program.  The 
description  of  control  characteristics  given  in  the  design 
report  was  used  to  make  a  qualitative  analysis  of  program 
operation. 

One  specific  condition  of  flight  was  selected  as  a  basis 
for  ccaparison  of  control  responses  to  various  combinations 
of  input  signals.  This  "base  condition”  of  flight  was 
defined  such  that  inputs  from  motion  sensors  could  be 
superimposed  on  base  condition  inputs  to  enable  an 
investigation  cf  the  effects  of  each  parameter.  The  base 
condition  for  program  testing  was  selected  to  model  an 
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TIME  (seconds) 

Fi9ar«  3.1.  B«ot  Ccaditioa  stabilator  Raaponsa 

aircraft  oparating  at  20,CC0  faet  and  250  tnots.  This  fixsd 

the  values  of  static  and  dyr.aBic  pressures  for  each 

siaulatlor.  and  peraltted  sanual  verification  of  the  gains 

produced  b;  the  air  date  schedules.  Additionally,  a  step 

function  representing  six  pounds  cf  force  in  the  aft 

> 

(positive)  direction  on  the  control  stick  vas  applied  at 
tine  0.0  seconds.  All  aoticn  sensor  inputs  were  held  at  zero 
so  that  their  effects  could  be  individually  studied.  The 
initial  deflections  of  the  stabilator,  leading  edge  flaps 
and  trailing  edge  flaps  were  set  to  zero  degrees  and  all 
initial  conditions  for  filters  were  set  to  zero. 


^•i  I'.oo  a'.a  t'.n  «'.a  s'.n  t'.oo  t'.qb  •’.«  *'.ao  ib.ot 
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Figtxre  3.2.  Bas«  Condition  Leading  Sdge  Flap  Besponse 

Figure  3.1  depicts  the  aoveeent  of  the  stabilator  in 
response  to  the  base  flight  conditions.  A  transient 
oscillation  is  produced  in  the  first  second  which  results 
oriaarily  froa  the  second  order  filter  in  the  pilot  stick 
dynaaics  path.  After  this  oscillation  has  decayed  the 
stabilator  continues  to  deflect  at  a  nearly  constant  rate, 
since  feedback  inputs  are  supressed.  The  stabilator  reaches 
the  liaiting  value  of  24  degrees  trailing  sdge  up  after  9.6 


seconds 
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Figore  3.3.  leading  Edge  Flap  Response  to  AOA  Feedback 


Figure  3.2  is  a  plot  of  leading  edge  flap  deflection 
(LEFIA?)  versus  "iae  for  the  bass  flight  condition.  A  steady 
state  flap  deflection  of  4.8  degrees  is  achieved  after  1.5 
seconds.  Since  lEFLAP  is  a  function  only  of  angle  of  attack, 
static  pressure  and  dynaaic  pressure  it  was  desired  to 
observe  the  variation  in  LEFLAP  with  angle  of  attack.  A  ramp 
type  increase  in  angle  of  attack  sensor  input  (AA1)  was 
superimposed  upon  the  base  flight  condition  beginning  at 
tine  1.4  seconds. The  resulting  schedule  of  leading  edge  flap 
deflection  is  shown  in  Figure  3.3  and  is  consistent  with  the 
functional  description  given  in  the  design  report. 
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Figure  3.4.  Base  Condition  Trailing  Edge  Flap  Response 

Trailing  edge  flap  deflection  was  modeled  similarly  in 
Figure  3.4  for  the  base  flight  condition  and  in  Figure  3,5 
for  the  condition  in  which  AA1  is  a  ramp  function  starting 
at  tine  1.4  seconds. 

To  determine  the  effect  of  forward  pressure  on  the  pilot 
stick  a  step  input  of  -12, C  pounds  was  superimposed  upon  the 
base  flight  condition  at  time  4,0  seconds,  which  simulated 
an  instant anecus  reversal  of  stick  force.  Figure  3.6  shows 
that  the  direction  of  statilator  deflection  changed  abruptly 


upon  introduction  of  the  new  stick  force. 


and  that  the 
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Figure  3.5.  Trailing  Edge  Flap  Besponse  to  kOk  Feedback 

magnitude  cf  the  steady  state  rate  of  stabilator  deflection 
was  approiiaately  unchanged.  Pilot  stick  force  inputs  of 
larger  magnitude  were  simulated,  but  the  results  are  not 
shown.  In  these  cases  the  rates  cf  stabilator  deflection 
increased  while  initial  oscillatory  behavior  exhibited 
characteristics  similar  tc  base  condition  rasponse. 

The  influence  of  angle  of  attack  feedback  on  stabilator 
oosition  is  displayed  in  Figure  3.7.  As  described  in  the 
design  report,  angles  of  attack  in  excess  of  22  degrees  will 
generate  a  proportional  stabilator  command  causing  the  nose 
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Figure  3.6.  Stabilator  Besponse  to  stick  Force  Reversal 

of  the  aircraft  to  pitch  down.  In  this  simulation  the  angle 
of  attack  input  function  AA1  was  set  egual  to  22  degrees 
plus  a  ramp  type  increase  of  one  degree  per  second  starting 
a-  time  tt.O  seconds.  The  resulting  stabilator  deflection  was 
identical  to  that  for  the  base  condition  until  time  4.0 
seconds  due  to  the  action  of  the  -22  degree  bias  in  the  AO& 
feedback  path.  The  nose  dcwn  pitching  effect  of  AOA  feedback 
is  apparent  at  all  later  times.  In  response  to  the  base 
condition  the  stabilator  had  reached  its  maximum  limit  of  24 
degrees  trailing  edge  down  at  time  9.5  seconds.  The 


40 


application  of  AOA  feedback  restricted  stabilaxcr  deflection 


to  20.9  degrees  at  time  9,5  seconds. 


rate  beginning  at  time  4.0  seconds.  The  change  in  direction 
of  stabilator  Bovement  is  much  more  rapid  than  that  produced 
by  angle  of  attack  or  norial  acceleration  damping. 


Figure  3.9.  Pitch  Rate  Damping 


The  thesis  computer  program  used  a  significant  portion 
of  allowable  CSHE  program  size.  Table  IV,  which  is  extracted 
from  the  "translation  table"  section  of  CSHP  output,  lists 
the  areas  of  the  program  which  most  closely  approached  the 


size  limits  of  CSHP 


lABLE  IV 


Program  Size  Restrictions 


Parameter 

HACBO  and  statement  outputs 
Statement  input  vorA  area 
Parameters-funct ion  generators 
History  and  memory  block  names 
BACRO  statement  storace 
SORT  sections 

Naximus  statements  in  section 


Current 

Program 


C5MP 

Maximum 

600 

1900 

400 

50 

125 

20 

600 


m  m 


The  computer  program  developed  in  this  thesis  simulates 
ihe  performance  cf  the  longitudinal  flight  control  system  of 
the  F/A-10  aircraft  by  generating  control  surface 
deflections  and  cross  axis  electrical  signals.  The  responses 
to  AOA,  normal  acceleration  and  pitch  rate  feedback 
correspond  xc  the  descriptions  given  b/  the  aircraft 
manufacturer. 

The  use  of  CSHP  simplified  the  task  of  system  modeling 
by  provision  cf  pre-programmed  functional  blocks  and  a 
flexible  format  for  outputs.  The  CSMP  orogram  size 
restrictions  do  not  appear  to  be  a  factor  which  would 
prohibit  the  addition  of  the  lateral  and  directional  flight 
control  systems  to  this  simulation.  It  is  recommended  that 


future  expansion  cf  this  prograa  be  done  using  techniques 
conserve  progran  size.  FORTRAN  subroutines  and  computer 
library  functions  should  be  utilized  due  to  the  limited 
MACRO  and  MOSORT  capabilities  of  the  CSHP  translator. 


tt5 


PITCH  RATE  gyro  PATH 


FORWARD  INTEGRATOR 
PATH 


_ ♦ 

^ITIVC  LlMlTERl 


FI4 


rx5 

_  _ 

OICITM.  FJLTCR 

R9 


FH 


V 


MIAIU  LIMITCR 


STABILATOR  PATH 


LEADING  EDGE 
FLAP  PATH 


TRAILING  EDGE  FLAP  PATH 


*tgBHDIX  E 


FLIGHT  COHTBOL  COMPUTER  PROGRAM 


*  **  ** 

*  *♦  LONGITUDINAL  FLIGHT  CONTROL  LAN  SIMULATION  ♦♦ 

*  **  ** 
*  *****************************  **********************  ***** 
* 

*  ************************** ************************* 

*  ♦♦♦♦♦INTEGRATING.LAG.AND  NOTCH  DIGITAL  FILTERS***** 

*  *************************************************** 

MACRO  FOUT-ZINT  (FIN  ,  K A,  KB  ,IMP, POUTZ  1) 

PROCEDURAL 

IFdMP.NE.  1.0)  GO  TO  10 
IF(KEEP.NE.  1.0)  GO  TO  10 
IF(TIME.EC.O.O)  GO  TO  10 
IODT*KA*PIN  ♦  KB*F00TZ1 
FOUTZ1«POOT 
10  CONTINUE 
ENDHAC 

MACRO  POOT*ZLAG(FIN,KA,KB,KC,IMP  ,PINZ1  ,FODTZ1) 

PROCEDORAI 

IF(IMP.NE.  1.0)  GO  TO  10 
IF  KEEP.  NE.  1.0)  GO  TO  10 
II(TIME.EQ.O.O)  GO  TO  10 
FOOT»KA*FIN  -  KB*PINZ1  ♦  KC*F00TZ1 
fOUTZ 1»POOT 
FINZIaFIN 
10  CONTINUE 
ENDHAC 

MACRO  FOUT*2NOTCH(FIN,KA^KB^KC,KD,KE,IMP^FINZ1,.  .. 
FIHZ2,FOUTZ1  ,FOOT22) 

PROCEDURAL 

IF(IMP.NE.  1.0)  GO  TO  10 
2F(KEEP.NE.  1.0)  GO  TO  10 
IFjTIME.EC.0.0)  GC  TO  10 

fODTaKA*FlN+KB*PINZ1+KC*FINZ2-KD*FOUTZ1-KS*FOOTZ2 
FOOTZ2  a  PC0TZ1 
IOOTZ1  a  FOOT 
PINZ2  a  FINZ1 
IINZ1  a  FIN 
10  CONTINUE 
ENDHAC 
* 

*  ***************************** 

*  ♦♦***FREOOENCI  AVERAGERS***** 

*  ***************************** 

MACRO  Z40aAV2040  (Z20,IHP) 

PROCEDURAL 

IF  (KEEP. NE. 1.0)  GO  TO  20 

IF  (TIME. HE. 0.0)  GO  TO  5 

Z20Z1=Z20 

Z40Z  1aZ20 

DEL«0.0 

GO  TO  10 

5  IF  (IMP. EQ. 1.0)  GO  TC  10 

Z40aZ40ZU0EL 
GO  TC  15 
10  Z40*Z20 

DEL*  iZ20-Z20Z1)  /2. 0 
Z20Z  1*Z2C 
Z40ZlaZ40 
15  CONTINUE 

20  CONTINUE 
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4 

\ 

1 


AC 

0  Z80«AV4080(Z40,IHP) 

fBOCEDOBAl 

IP  (KEEP. NE. 1.0)  GO  10  2( 

IP  JTIBE.HE.O.O)  GO  TO  5 

Z40Z1«Z40 

Z80Z1-Z40 

DEL»0.0 

GO  TC  10 

IP  (IMP.EQ.  1. 0)  GO  TO  10 

Z80«Z80Z1>CEL 

GO  TO  IS 

Z80«Z40 

DBL*iZ40-Z40Z1)  /2.0 

Z43z1>Z40 

Z80Z  1«Z80 


15  COHTIHOE 
20  CONTIHOE 
END  N iC 

HACBC  Z80-AV2Q8C  (Z20,IHP) 

EBOCECOBAl 

IP  (REEP.NE.1. 0}  GO  TO  20 
IP  jlIME.ME.0.0)  GO  TO  5 
Z20Z  1«Z20 
Z80Z1>Z20 
GO  TO  10 

5  IP  (IBP. EQ. 1.0)  GO  TC  10 

Z80«Z80ZUDEL2 
GO  TC  15 
10  Z80SZ20 

DEL2«(Z2C-Z20Z1)/4.0 
15  COHTIHUE 

Z20Z1»Z20 
Z80Z  1«Z80 
20  COHTINOE 

EBON  AC 
* 

♦  ♦♦♦♦♦INITIAL  CONDITIONS  AND  CONSTANTS  FOR  BASE^^^^^ 

♦  ♦♦♦♦♦PLIGHT  CONDITICN  OP  2  0,000  FT  AND  250  KTS^^^^^ 

*  *«««*«**«*4>4i4i*««4i4>«***4i***  « A*  ««« «««««*  *«*««« 

IN  ITIAL 

CONSTANT  PS»972. 49.  QC*1 12.73,  RAT  E20*0. 05.R  ATE40»0.  025,  . . . 

PS 3IC1*0.0,PS3IC 2=0. 0, PS4IC1«0. 0  ,PR2IC  1=0.0, PR2IC2*0.0  ,.  .  . 
AA2IC2=0. 0  .NZ2IC1«0 . 0 ,BZ2IC2«0. 0 ,NZ4Z1=0. 0 ,NZ5Z 1*0 .0, . . . 
S5.^o«Oiii^^/»f|6z^!?^A2A?P^S^*9*o»8P!’z2=o.6,HPi2zi=o.o,... 
P8A  =  0.69084-F8E»-0.9  9068.  E8C»0. 6  6312,  P8D*-0l  99068,  ... 

DS  1*0.0,DLE1=0.0,DTE1»0.  b,LE2Z1«0.0,IE2Z1=0,  0,LEDIC1*0 . 0  , . . 
LEDIC2=0.O-TEDICi«O.0.TEDIC2«O.O,PR3Z1*O.O,PR4Z1=O.O,.  .  . 
RATE80*0.0125,AA2IC1«0.  0,  EK9«-1 .  1543,  HP12 Z2=0 . 0  ,P8  S»0.  3539 

DTNAHIC 

* 

SORT 

♦ 

PS1«6.0^STEP(0.0) 

PR1«0.0 

AA1>22.0>BAIIP  (4.0) 

NZ1«0.0 

RR1-0.0 

TR1»0.0 

* 

* 

*  *4i«««**««*******«««4«*«««*«««***4i4i4i* 

♦  ♦•♦♦♦PONCTIONS  REQOIRING  NOSOHT^^^^^ 

«  *4141  «*«******«««*«%« «««««**  •4i««4i4i*4i4i« 
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r.  i.".  •r.  V. 


M 

I 

M 


¥ 

m 


J 


L-V- 


L^» 


* 

« 

* 

m 


* 

4> 

* 

« 


««****  «««****«**4i««  44 

♦♦♦♦♦FONCTION  12***** 

4t «* 4  *4  44 4 « 4 44t  4  *  44<* 41 4  4 

F12?1»Il**2*9.6  25-.025*SI+  1.0 
Fl2T2=PS*7.969E-44  0.eu 
F12HAX»lIHITn  .0.8.  0,F12T21 
F12T3*LIHIT  (1 . 0  ,F  1 2 MAX,  FI 2T1) 
F12T5«LIMIT(.5,  1-35. HI) 
F12T6aF12T5«(0.00952*PS+4.  04) 
F12T4»LIBIT  (1.0.8.0,F12T6) 
NCSORT 

IFi.HOT.BI.GT.0.5)  60  TO  5 
F12  *  P12T4 
60  TO  6 
CONTINOE 
F12  »  F12T3 
COMTINOE 
SORT 

444444444444444444444 

4444  4P0NCTIOII  40***** 

444444444444444444444 


♦  (-PS*0.  00396-  1.  18) 


F40T  1*3- 25-3.  0*LIMIT  (0.75,  0.8 5. HI) 
PaOT2=0.65625-.0013  125*LIMIT(56o.O,1800.0,PS) 
P40T3*-0.26177+ (9.635E-4)  *LIMIT  (500 .  O.  1800. 0, PS) 
PIQ*III!1T{0.0, 1800.  C.PS) 

P40T4alIBIT  (0.0,105.0,00-3  35.0) 

F40T5*  (F40T4*  (-PIQ*  (1.17428E-6)  ♦  (1.5238E-3))  )  ... 
♦  0.475-]6.5E-4)4piq 


F40T6»F40T34(F40T2*LIHIT{0.75,0.85,RI)  ) 
P40T7»P4flT44j5.6746E-44l.984lB-7*PIQ)  .  .. 
♦  JO.  16923-3. 84615E-5*PI9) 

(0.  16923-  3.  86415E-5*Pi5 

—3  -9.2 


F40T8=(  _ _  _  _ 

♦  (F40T4*(  1.67247b 


'152E-74PI0) ) 


20 

30 


40 

50 


60 


70 

80 

4 

4 

4 

4 

4 

4 

4 


NCSORT 

IF  (. NOT. RI.GE.  0.75)  60  TO  20 
F40T9*F40T8*F40T1 
GO  TO  30 
CONTINOE 

F40T9*F40T8 

CONTINOE 

IF  (.NOT. (OC.GE. 440. C. AND. BI.6E. 0.75) )  60  TO  40 
F40T10»F40T6 
GO  TO  50 
CONTINOE 

F40T10»F40T7 

CONTINOE 

IF  (.NOT  .PIQ.6T  .980.0  60  TO  60 
F40»F40T9 
GO  TO  80 

CONTINOE  , 

IF  (. NOT  .PIQ.LE. 500.0  60  TO  70 
Fa0»F40T5 
GO  TO  80 
CONTINOE 

F4C«F40T10 

CONTINOE 

*44  4  4  4  4  4  44  444  4  4  44  4  4  4  444  444 
44444SOHTED  FONCriOHS4*444 
444  44444  44  44**444*4444*4*4 

SCBT 

444 4 44 44444*4 4 4 *4 *44 4*4 44 
44444P0KCTICNS  22,2  3*4*4* 

444  444  44  444*44  444  4  4  4444  44 

P22*0.01674(-80  0.0  +  lIl!IT(8  00.  0,  900.0, QC) ) 
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*»»  ****  ♦***  *  ****  *  ******  **** 


P23«  3.  14  35  -0.  14  29*LIB1T  {15 . 0, 21 . 998 ,  &LPH&T) 


♦♦♦♦♦PONCTIONS  25,27.28,29  0,29L,32A,37***** 

**«*««*«*«***«*««4i«*A*«4i*«  «****4ii^  ' 


F25»47. 636-0. 05106*LIHIT  (6  00. 0,835.0,00 


P27»1.328*  (ALPHAT+7.6584-l7.86*LIMIT  (0.44 ,0.  63,RI)  ) 
P28»44.5£1-0.0  4  05  8«IINIT(2  60.  0,950.0  ,QC) 


P290«87.3825-76.25*1ISIT(0.7,  1. 146,RI) 


P291«0.0 


? 


KF-lIHII  1200. 0,2000. 0,QC) 
32Aa100.0/QKP 


P37a2.5-0.  5*LIMIT(3.C,5.0,  NZA) 


*****«««  4*  ♦^i************** 

♦♦♦♦♦PONCTIONS  68,107***»* 

P68«-0, 002977*  (-480. C*LIHI  T  (250 . 0,48 0.  0,QC) ) 


P107T1*  1  (-5.714E-71  *ABS(QC-75C 
P107«LieiT(0.  0,  10000. 0,P107T1) 


C-750.0))  ♦  (8.4E-4) 


«*««**«« 4********** 4 4^****  * 

*****I8P0LSE  POTICTICliS***** 

**«««««*4*«4>*«*4***44****** 

iaP20«IHF0LS(0. 0,BA7E20) 
iaP40»IEP01S(0.0,RAlP40) 

iap80»iepais(0. 0, baheoo) 


*44  *4i*««4*««**  ******** 
44444PK  FOBCriOHS***** 
44******4***4*****4*4* 

PK1 1«1.65*F22 
PK12«0.5654 
PK16A«S.5*F32A 
PK17A*P37*F32A 
PK19A*f  12*F32A*3.5 
PK21A*P12*F32A*F68 
PK22A«Pl4«F32A 
PK16«ZHCID  (IHP20,PK16A) 
PK17«ZBC1D  (iaP20,PK17A) 
PK19»ZHCI0  |iaP20,  PK19A) 
PK21»ZBOLD(IBP20,PK21A) 
PK22»ZBCI0  (iaP20,  PK22A) 


44**4*******4* ************ ****** 

*****PI10T  STICK  IHIOT  PATB***** 

44**4**4****4*************  *4**4* 

PS2-CEAESP  (-2. 0,2. 0 ,FS1 ) 

PS3*CaPXPl  1PS3IC1,PS3IC2,0 . 14, 27.  3,  PS2*  106 .47) 
PS4«BEAlPliPS4ICl,7.S365E-3,PS3) 

PS5-ZHC1D  (IBP80,PS4) 

PS6-PS5 

PS7«PS6* J7.0*  (0.2*AES(PS6) 

PS8«liaiT(-25.0,50.  0,PS7) 

PV3«FS8 

P73A»Z BOLD  (IRP20,PV3) 

PS9«F?3*F32A 


)) 


****************************** 
*****PITCB  bate  gtbc  path***** 
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m 


4>4i  4>  4  *4i*4>**4i4>4i4>4>4t4i4i  *4**41*4141  **** 

PR2»CBPXEl  (PH2IC1,PB2IC2,  0.89, 78.5,  P  81*6162.  25) 
PB3«ZBOID(II1P80  ,PR2) 

P2A»1.0*EK  11*(1.5-PK12) 

P2B»i1.0+PK11)  *  (1.0-EK12) 

P2C»1.0-PK12 

PR4>ZL&G  (PB3.P2A,  P2E.P2C,inP80,PR3Z1,P8az  1) 


************************************* 

*****j^II61E  OP  ATTACK  SENSOR  PATH***** 
************************************* 

AA2>CHPXPL  (U2IC1.  A  A2IC2,  0.74, 209.0,  AA1«4  368  1.0) 
AA3SZH01D  (IBP40,AA2) 

ALPBAT*0.S9*AA3-H.9 
ALPHAS»AIPHAT 
AA6-ALPBAT-22.0 
AA7«1I(]IT  (0.0, 10000. 0,AA6) 

P?1*PK17*AA7 
PT1A«ZHC1D  (IHP20,  PV1) 

*********************************** 

*****M0BBAI  ACCELERCHETER  PATH***** 
*********************************** 

NZ2«CaPZPL(IlZ2lC1,HZ2IC2,  0.89, 200.0,  HZ  1*4000 0.0) 
NZ3«ZB01D  (IBP40,HZ2) 

NZA«NZ3 

HZ4»HZ3-  (BB1**2)*  (6.8529E-6) 
P5A«(1.0*PK9*(1.0-PK10)  ) 

P5B*  j1.0+FK9)  *(1.0-EK16) 

P5C»1.0-PK10 

MZ5»ZLAG  JNZ4.PSA.PSE,P5C,IHP40,HZ4Z1  ,NZ5Z1) 

PV2*lIHITi(-10.0,lO.O,HZ5f 

»Z7*P?2*EK16 

PV2A«ZBOlO(inP20,PY2) 

********************************* 

*****fOBfiABO  INTEGRATOR  PA  IH***** 
********************************* 

EI2>HP3 

PI3»IIBIT  (-10  000.  0.0.0,  FI2) 

PI4*LIHIT(0.0,  10000. 0,PI2) 

FI5»  (FI3*F23)  ♦PW 
FI6*ZIHT(PI5,0.  05,1  .C,IHP2  0,PI6Z1) 


******************* 

*****HAxN  path***** 
******************* 

aP1*PV3A*PK22 
HP2«  (P?4A*PK21)  -HPI 
HP3*  (PK19*E?2Ai  ♦PTI  A+HP2 
HP4*LIHiT  (-50. 0  ,  25.  C,FI6) 
HP5A«A72040  (HP4  ,IHP20) 
HP5»ZHOIO  (IHP40  ,HP5A) 
HP6*HZ7*KE5 
HP7*aP6*EV1 

HP8A-HP7*  (  (RR1*TR1)  ♦F107) 
nP8*ZHOIDjlfIP40,HP8M 
aP9»AV408d  (aP8,IHP4  0) 


HP8A-HP7*  (  (RR1*TR1)  ♦F107) 
nP8*ZHOIDjlfIP40,HP8M 
aP9»AV408d  (aP8,IHP4  0) 

HP10«BP9*PB5 

HP11»LIHIT  (-25.0,25.C,HP10) 

aP12»ZNCTCH  (HPI  1.P8A,P8B,P8C,P8D,P8E,iaP80,.  .. 
HPIIZI^PI  1Z2,aP12Z1,  HP12Z2) 


ar  I 

HP13-HP12-PS9 


************************* 
*****STAEIlArOR  PATH***** 
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*  4i*«**«*«4««4t*****««**«4>*« 

ST2«lIf!IT(-24.0,10. 5,HP131 
ST3*ST2-DS1 

ST4«LIHIT(-24,0,10.5,ST3) 
ST5*C1*TZBJ0.012  5-ST4) 

^  STiDfF«CBPXPL  (TE6iC1,TEDIC2, 0.7, 40.0 

«  *«4>**««*  4* ««*««««««*««««««  **«4>*« 

*  ♦♦♦♦♦LEADING  EDGE  FLAP  PATH^*^^^ 

«  «*«**4>**««*Mi*«***«*««*««*****4i« 

LE2«ZINT  (ALPHAS. 0.0625, 0.9 375, IHPSO, 
LE2T«ZHCI0  (IHP20,LE2) 

LE2T1A=II2T^F27 

LE2T1>ZBOLD(IHP20.LE2T1A) 

LE2T2alIHIT  (0. 0  .F28  .LE2T1) 
LE3»LIHIT(C.0,F29O,£e2T2) 

4t  •^^^iik*^*****  *4i4i*4i4>4t4t4>* 

«  ***«*18  DEG/SBC  RATE  LIMIT, 

«  ^t^^i^^^^t****’******  ********** 

HO  SORT 

IF  (TIME.  »E. 0.0)  GO  TC  85 
LE4Z1SLE3 
85  COHTIMOE 

DEL  1*IE3-LB4Z1 
DELIie=lIMIT( -0.9, 0.9,  DELI) 
LE4A*LE4Z1<»DELLIH 
IF  (IMP20.HE.  1.0)  GC  TO  90 
LE4Z  1SLE4A 
90  COHTIMOE 

LE4>ZHOLDiIHP20.LE4A) 

*  *******ii********** ********** 

SORT 

LE5«AV2080  (LE4,IHP20) 

DLE2P1»LE5*3.0 

DLE2-DEADSP  <-0LE2Pl  ,CLE2P1  ,  DLE1) 
DLE3*DLE1+DLE2 

*  H0TE:P0SITI7E  LEADIHG  EDGE  FLAP  PATH 
LE6»IE5<»DLE3 

LB7sIIMIT  (-3.  0,33.0. IE6) 

*  ****:************************ 

*  ****0^8  DEG/SBC  RATE  LIMIT, 

HO SORT 

IF(TIHE.NE.O.O)  GO  TO  110 
LE8Z1-IE7 
110  COHTIMOE 

DEL1«IE7-IE8Z1 

0ELLIM>LIHIT(-0. 225,0. 225,  DELI) 
LE8A«IE8Z14DELLIH 
IF  (IHP80.HE.  1. 0)  GO  TO  12  0 
LE821«LE8A 
120  COHTIMOE 


,ST5*1600.0) 


LE2Z1) 


********** 

20 

********** 


*********** 


OMLI 


********** 

80 

********** 


SORT 


LE8«ZRCIDJIRP80  .LE8  A) 

«««*«*«*4>«**4 


*«««****  ****«4(4c«***4c*4i4(4>*««4t« 


LE9«LIHIl(-2.  0,33.0.1E8) 

LEI  0«CHTZB  <0.  0 1  25,  LE9) 

LEFLAP«CHPXPL (LEDIC 1  ,LE0IC  2, 1.4,20.0 


,LE10^400.0) 


4i*4i*««««««*«*««««4i«*******  *** 

♦♦♦♦♦HOSCHT  FOMCTIOH  24  ♦♦♦♦♦ 
***************************** 

P24 I 1-22.538-20. 51 ♦LIMIT (0.27,0. 66, R 
F24L2-32.76-36.0^LIHIT(0.6  6,0.91,RI) 
F24T1*LIHIT(0.  0, 10000.0, ALPHAT) 
P24T2-AIEHAT-  (1  4.  8769-7.69  23^LIMIT  (0 
F24T  3*- 2.  DELIMIT  (0.0, 10000. 0,P24T2) 
F24T4-  1.4^  (P24T  1-fF24T3) 

HOSORT 


.27,0. 91, RI)  ) 


IP 


0.66)  GO  TO  10 


F2ai»P24l2 
GO  TO  15 
COHTmOE 

F24L»F2aL1 

COMTINOE 

F24sLIMIT(0.0,F24L,F2aT4) 


PAGE  XTPLOT 
OOTPOT 
END 

STOP 

ENOJOB 

&  ■ 

4t  «*««***««*«««**«**«****«« 

*  ♦♦♦♦♦TBilLIHG  EDGE  FLAP  PATH***** 

*  4t«**«««*j»«4iMi**»«**«***«4i«  «*«*«*« 

TE2«ZINT  (ALPHAS. 0.03 125,0.  096875  , IBP  80  ,TE2Z1) 
TE2T»ZHCLD  (IHP20,TE2) 

TE2T1«TE2T*P24 

TE3«LIMIT(0.0,F25,TE2T1) 

*  NOTE:POSITIVE  TBAILIUG  EDGE  FLAP  PATH  ONLY 

TE4»TE3*DTE1 

rE5»lIi!IT  (-8.0.45.0.TE4) 

TE6«QNTZHi0.0l25.TE5) 

TEPLAE=CBEXPL  (T  EDICI  ,TEDIC2, 1 . 4, 20.  0  ,TE6*40.  0) 
TEBHINAL 
HETHCD  BKSFX 

TIHEB  PIHTIH»10. 0,OOTDEL«0. 0125,  PBDEL»0. 1250, DBLT=0. 01  25 
PBIHT  STAEEF,LEFLAP,TEFLAP 
LABEL  STADEF  IN  DE6BEES 
LABEL  TINE  IN  SECONDS 


TEDIC2,1.4,20.0,TE6*40.  0) 


AIR  CATA 


SCHEDULES 

LABEI  FUNCTION  12 

* 

INITIAL 

PASAHETER  ?S-  (2000.0,100  0.0,500.0,230.0) 

* 

DTNAHIC 

QORT 

RIsRAHP  (0.0) 

P12T1»RI**2*9.625-.  025*HI*  1.0 
P12T2apS*7.969E-4+0.e4 
F12HAX=IIHIT(1 .0,8. 0.F12T2) 

F12T3=LIMIT  (1 . 0  ,P  12MAX,  F12T1) 

F12T5«LIMIT(.5,  1.35 -HI) 

F12T6=P12T5*fO.  00  952*PS+4.  04)  ♦(-PS* 3.  00396-  1.18) 
F12T4*LimT  (1.0,8.0,F12T6) 

MCSORT 

IFj.NOT.HI.GT.O .5)  GO  TO  5 
F12  «  F12T4 
GO  TO  6 

5  CONTINUE 
F12  «  F12T3 

6  CONTINUE 

« 

TERMINAL 

TIMER  PINTIH  «2.0,  CUTDEL-0.05,  PRDEL=0.05 
PRINT  F12T1,F12T2,F12T3,P1  2T4  ,F1 2T5  , FI  2T6  ,P1  2MAX  ,  F1 2 
PAGE  XTPICT 
OUTPUT  HI,  F12 
END 
STOP 
ENDJCB 


LABEI  FUNCTION  24 

* 

INITIAL 

^  PARAMETER  HI* (0.27, 0.66,0.  91) 

DYNAMIC 

« 

SORT 

ALPHA*RAHP  (0.0) 

P24 1 1*22. 538-20.5 1  ♦LIMIT  (0.27,0. 66,  R  I) 

F24L2«32. 76-36. 0*LIMITi0.66, 0.91  ,RI) 

F24T1*LieiT(0. 0,10000. d, ALPHA) 

F24T2*ALPHA-( 14 .8769-7, 692  3*LIMITiO . 27, 0. 91, Rl) ) 
P24T3*-2.0^LiHIT(0. 3,10000 .0,F24T2) 

P24T4-1.4*  (F24T1+F24T3) 

* 

NOSOBT 

IF  J.MCT.HI.GT. 0.66)  GO  TO  10 
F^4L*F24L2 
GO  TO  15 
10  CONTINUE 

F24I»F24L1 
15  CONTINUE 

^  F24*IIHIT  (0.0,F24L,F24T4) 

TERMINAL 

TIMER  FINTin*30.0,O0TDBL*0 .5.PRDEL«0.5 

PRINT  F24L1,F24L2,F24T1,P2  4T2,F24T3,F24T4,F24L,F24 
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END 

STOP 

ENDJC6 


LABEL  FUNCTIOH  29 
* 

INITIAL 

PARAHETEB  PS«  (1250.0) 

* 

OYNAHIC 

BI>BAHP  iO.O) 

P290-87.  3825-76. 25*LIBIT  (0.7,  1.  146,81) 
P29L»0.0 

* 

TERHINAL 

TIHEH  PINTIH*1.5.  PRDEL*0.05,  O0TDEL=0.05 
PRINT  RI.P290,P29L 
PAGE  XYPtOT 
OOTPOT  BI,P290,P29L 
END 
STOP 
ENOJOE 


LABEL  F0NCTIC8  40 

* 

INITIAL 

PARABETER  QC*  (1700.  C,.440.  0 ,335 .0) 

* 

OYNAHIC 

SORT 

PS«  0.0001  ♦  RAHP  (0.0) 

F40?1»l?2S-3.  0*LIMIT{0.75,  0.85.RI) 
P40T2»0.65625-.0013  125*LIMIT  (500.0. 1  800.0. PS) 
P40T3»-0.26177+ /9.635E-4)  *LIMIT(500.  0,  1800.0, PS) 
PIQ»LIHIT]0.0, 1 800, O.PS) 

F40T4»LIHIT  (0.0, 105.0, QC _ 

P40T5»(P40T4*  (-PIQ*  (1.T742  8E-6)  ♦  (1. 5  238E-3) )  ) ... 
♦  0,475-j6  .5fi-4)*PIQ 

P40T6«P40T3+(F40T2*LiaiT(0.75,0,85iRI) ) 

F40T7>  “  '  '  ■ 


LIBIT  (0.0,105.0,QC-3  35.0) 

(PaOT4*  (-PIQ*  (1.T7428E-6) 

♦  0,475-j6.5fi-4)*PIQ 
-P40T3+(F4  0T2*LiaiT(0.75,0,w^,«*,  , 
»P4QT4*j5,6746E-4>1.9841B-7*PIQ)  .  .. 

♦  (0.  169  23-3. 84615E-5  ♦PIQ) 

«  (0.  16923-  3.  86415E-5*PIQ)  .  . , 

-------  - -3-9. 2r - 


P40T8«(  _ 

♦  (P40T4^  (  1.67247E 


!9152E-7^PIQ) ) 


20 

30 

40 

50 


NCSORT 

IP  (. NOT. HI.GE.  0.75)  GO  TO  20 
P40T9»P40T8*P40T1 
GO  TO  30 

CONTIHOB  _ 

P40T9»P40T8 

CONTINUE  ^  ^ 

IP  (.NOT.  (QC.GE.440.C.AND.  RI.GE.0.75)) 
P40T10*F40T6 
GO  TO  50 
CONTINUE 

P40T10-P40T7 

CONTINUE 

IP  (.NOT. PIQ.  GT. 980.0)  GO  TO  60 
P40«P40T9 


GO  TO  40 
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GO  TO  80 

60  CONTINOZ  ,  ^  > 

IP  (. SOT. PIQ.LE  .500.0  GO  TO  70 
F40»FUOT5 
GO  TO  80 

70  CONTINOE 

FaO=FUOTlO 
80  COSTINOE 

* 

TERM  INAL 

TIMER  FINTIMs2000.0.PRDEL»20.  0,OOTDEL=20. 0 
PRINT  PS,QC,RI,F40T7,FU0T9  ,F40T10,F40 
PAGE  XYPIOT 
OOTPOT  PS,P40 
END 
STOP 
ENDJCB 


iBEL  GENERATED  FONCTIOH 

THIS  FHOGRAM  GENERATES  QC  FDNCTIONS  22, 25,  28,  32A,  AND  68 
IITIAL 

PARAMETER  PS»1  000.0 
NAHIC 

QCsRABP  (0.0) 

P22-0.0167*  (-800.0  +  LlMIT(800.0,  900.0,QC)  ) 

F25»47. 636-0.  05 106*LIM1T(6  00.0, 835. 0  ,QC) 

P28=44.551-0.  04  058*LIM1T (2 60.0,950.0,00)  ' 

QKFalIMITj(200.  0,2000.0.QC) 

F32i»100.(J  /  QKF 

F68  =  -0.002977*(-480.0+LIMIT  (260. 0,480.0,00))’ 

IRHINAL 

TIMER  PINTIH-2500. 0.  O0TDEL»50.0,  PRDEL*50.0 
PRINT  GO  ,P22,F25,F28,QKF,F32A,P68 
PAGE  XIFIOT 
OOTFOT  CO,F22 
OOTFOT  QO,F25 
OOTFOT  CC,F28 
OOTPOT  QO,F32i 
OOTFOT  CC,F68 
END 
STOP 
IDJOE 


LABEL  GENERATED  FONCTION 
INITIAL 

PARAMETER  RI> (0 .44 , 0.63) 

* 

DTNAHIO 

ALPHJI»  SAHF  (0.0) 

MZA«AIPHA 

ALPRAT«IIPRA 

* 

F23-3.  1435-0.  1429*LIHIT  (15 .0,  21.998,  ALPHAT) 

* 
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F27a  1.328*  ( AL PH  A+7.  8584-17  .86*LIMIT  ( 0. 44, 0. 6  3, HI )  ) 
F37s2.5-0.5*LIMIT  (3 . 0  ,5.0  ,  NZA) 


TEBHINAL 

TIMEB  PINTItt=30.Q, 
PRINT  F*3,F27,P37 
PAGE  XYPIOT 
ODTPOT  ALPHAT,F23 
OUTPUT  AIPHA.F27 
OUTPUT  NZA,F37 
END 
STOP 
ENDJCB 


EBDEL=1.0,  OOTDEL»1.0 


AgPENDI X  D 


COMPOTES  PROGBAMS  FOR  SIGNAL  BLOCK  TESTING 


CEADSPACB  PONCTION  AND  ALIASING  FILTERS 


INITIAL 

CONSTANT  ES2C1*2.0 

♦ 

DYNAMIC 

PS1»  -5.0  ♦  RAIIP(O.O) 

PSPaDEADSP  (-PS2Cl-PS2C1,PS1) 

PS3-CMEXP1  (0.0, 0.5.0.  14.27.3  .PS2*106. 47) 

^  PS4aHEALPL  (0.0, 7.  9365E-3,PS3) 

TERMINAL 

TIMER  FINTIM=10.0,  OOTCEL=0. 25,  PRDEL=0.25 
LABEL  STICK  DEADSPACE  PLOT 
PRINT  ES1,ES2 

ODTPOT  TIME,  PS  1 , PS2 , FS3, PS4 
PAGE  XYPLCT 
END 
STOP 
ENDJCE 


FBEQOENCY  SVEEAGERS 

MACRO  Z80  =  AV2080  (Z20,IHP) 
FROCEDORAl 

IP  (TIME. NE. 0.0)  GO  TO  52 
Z20Z1*Z20 
Z80Z  13Z20 
GO  TO  54 

52  IF  (IMP. EQ.  1.0)  GO  TO  54 

Z803Z80ZUEEL2 
GO  TO  56 
54  Z80>Z20 

DEL2»{Z2C-Z20Zl)/4.  C 
56  CONTIHUE 

Z20Z  1»Z20 
Z80Z1=Z80 

ENDMAC 

INITIAL 

CONSTANT  HATE»0.05 
DYNAMIC 

Y1=  RAMP  (0.0) 

YIAaEXP  (Y1) 

Y2*  IMPDLS  (0.  O.RATE) 

Y3*  2HOIDJY2.Y1  A) 

Y4=  AV2080  (y3,  Y2j 


TERMINAL 

METHOD  RRSFX 

TIMER  FINTIB»1. 0,OOTCEL  = 
LABEL  20  TO  80  AVERAGER 
PRINT  Y2,Y3,Y4 
END 
STOP 
ENDJOE 

MACRC  Z40«AV2040  (Z20  ,IHP) 
PROCEDORAL 

IP  (TIME. NE. 0.0)  GO  TO  5 

Z20Z  1=Z20 

Z40Z1sZ20 

DEL=0.0 

GO  TC  10 

5  IP  (IHP.EQ.1. 0)  GO  TO  10 

Z40»Z40Z14-DEL 
GO  TO  15 
10  Z40»Z20 

DEL*  (Z20-Z20Z1)  /2.  0 

Z20Z1*Z20 
Z40Z  1»Z40 
15  CONTINUE 


0,OOTCEL=0.  0125,PRDEL*0.0125,  DELr=0.  012  5 
AVERAGER 


BA  TE) 


15  CONTINUE 

ENDMAC 

INITIAL 

CONSTANT  RATE*0.05 
DYNAMIC 

Y1=  HAMP(O.O) 

YlAaEXPJll) 

12*  IMPUlsrO.  O.RATE) 

13*  ZHOIE  (Y2,T1  A) 

Y4»AV2040(Y3,  Y2) 

TERMINAL 

METHOD  BRSFX 

TIMER  PINTIM*1. 0,OaTDEL=0.  025  ,PRDEL  =  0 . 025 » DBLT*0 . 01 25 
LABEL  20  TC  40  AVERAGER 
PRINT  Y1  .Y1A,Y2  ,Y3,Y4 
PAGE  XYPLOT 
OUTPUT  TIBE,I1A,I3,T4 
END 
STOP 
ENDJOE 


BATE  LIMITERS 
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DIGITAL  FILTEBS 


BACHC  POOT«ZIIIT  (FIN,  K A, KB  ,Il!P,FOOTZ1) 
FBOCEOOBAl 


IFflHP.NE.  1.0] 
IFfKEEP.NE.I.a 
IIjfTIHE.EQ.O.Ol 
FCOT-KA'i'PIK  ♦ 
FOUTZ 1*FOOT 
CONTI NOE 


GO  TO  10 
GO  TO  10 
GO  TO  10 
B*POOTZ1 


ENDHAC 

HACBC  FOaT>ZLAG(FIN,  KA,RB,KC,INP,FINZ1,POUTZ1) 

PBOCEDORAL 

IF(IHP.NE.  1.0)  GO  TO  10 
IFfKEIP.NE.  1.0)  GO  TO  10 
IF(TIJJE.EQ.O.O)  GO  TO  10 
IOUT»KA*FIN  -  KB*FINZ1  ♦  KC*F0DTZ1 
FOOTZI-POOT 
IIN21»FIN 
10  CCNTINOE 
ENDHAC 

HACBC  FOOT*ZNCTCH  (FI  N  .KA,  KE.KC.K  E,KE,IHP,  .  .. 
FINZ1,nNZ2,FOCIZ1,FO  0TZ2) 

PBOCEDOBAl 

IF(IHP.NE.  1.0)  GO  TO  10 
IF  KEEP. NE. 1.0)  GO  TO  10 
IFfTIHE.EQ.O.O)  GO  TO  10 

FOOT=KA*FIN+KB*FINZ1+KC*FIN2  2-KD*FOaTZ  1-KE*FO[IT22 
FOOTZ2  *  FOOTZ1 
FOOTZ1  »  FOOT 
FINZ2  *  FIN21 
FINZ1  «  PIS 
10  CONTINOE 
ENDHAC 
INITIAL 

PA  BAH  INTN».0125.INTD»1.0 ,AA1 1A1 «0. O.AAIIAsO.O 
PARAH  PK9*-1.  15<43,  PK  10-0.4647,  A  A 1 0BUO.  0  ,  AA 1 1 B  1«0 .0 
PARAH  NA». 69084. NB—.  99068-NC-.6631  2,ND=-.  99068, NE=.  35  396 
PABAH  AA1GC1-0.0,U10C2«0.0,AA11C1«0.0,AA1  1C2-0.0 
LA-(1.+PK9*(1.-FK10)  ) 

LB-(1.+PK9)*(1.-PK10f 

LC-1.-PK10 

DTNAHIC 

AA9A  >  STEPiO.O) 

IHP-IHPOLS  (0.0,0.01  25) 


IHP-IHPOLS  (0.0,0.01  25) 

AA10A  -  ZHCID  (IHP,AA9A) 

AA11A  -  ZINT(AA10A,INTN,INTD,inP.AAllAl) 
AA11E-ZLAG  (AAIOA.LA.LB,  Lc ,IHP,  AA 10B 1 ,  AA1 1 B  1) 
AA11C  -  ZNOTCH(AA10A,NA,NB,NC,ND  ,NE,inP,AA10( 
AA10C2,AA11C1,AA1  1C2) 

TERHINAL 
HETHCD  RKSFX 

TIHEB  FINTIH»1.5,OaTDEL-0.0125,DELT-0.0125 
PRINT  AA10A.AA11A.AA1 1B. AA11C 
OOTPOT  AA10A  (0.0,^5)  ,AA  1 1A  (0.0,  1 .5)  ,  . . . 
AA11Bb.0,1.5i  ,AA11C(0.0,  1.5) 

EN  D 
STOP 
ENDJOB 
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